Quantum Confinement Induced Magnetism in LaNi03-LaMn03 Superlattices 
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The emergence of magnetic reconstructions at the interfaces of oxide heterostructures are often ex- 
plained via subtle modifications in the electronic densities, exchange couplings, or strain. Here an ad- 
ditional possible route for induced magnetism is studied in the context of the (LaNi03) n /(LaMn03) n 
superlattices using a hybrid tight-binding model. In the baNiC>3 region, the induced magnetizations 
decouple from the intensity of charge leakage from Mn to Ni, but originate from the spin-filtered 
quantum confinement present in these nanostructures. In general, the induced magnetization is the 
largest for the (lll)-stacking and the weakest for the (OOl)-stacking superlattices, results compatible 
with the exchange bias effects reported by Gibert et al. Nat. Mater. 11, 195 (2012). 

PACS numbers: 75.70. Cn, 75.75.+a, 75.47.Lx 
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Introduction. The area of research that focuses on ox- 
ide heterostructures is attracting considerable attention 
because of its importance in the development of quantum 
devices based on correlated electronic systems Novel 
physical properties are expected to emerge from the elec- 
tronic reconstruction near the interfaces^— In particular, 
the interfacial magnetism can have properties different 
from those of bulk materials and in recent years several 
investigations revealed various magnetic reconstructions. 
Their origin can be mainly classified via mechanisms in- 
volving modifications in the (1) electronic densities, (2) 
exchange couplings, and/or (3) strain. For example, 
mainly due to strain and charge transfer, pure mangan- 
ite LaMn0 3 /SrMn0 3 (LMO/SMO) superlattices (SLs) 
display a variety of magnetic orders ^r— while due to 
modifications in the exchange coupling, distinct magnetic 
states emerge in the (001)- (011)- and (lll)-stacking of 
LaFe0 3 /LaCr0 3 SLs^~— Furthermore, the interfacial 
magnetic orders of Lai-^Sr^MnOs (LSMO) can be tuned 
by attaching ferroelectric layers (e.g. LSMO/BiFe03, 
LSMO/BaTi0 3 , and LSMO/PZT), that induce modifi- 
cations in the interfacial electronic density*^ - — 

Recently, an exchange bias effect was reported exper- 
imentally in (LaNi0 3 ) n /(LaMn0 3 )„ (LNO/LMO) SLs 
grown along the (111) axis of a pseudocubic structured 
In contrast, no exchange bias was observed in the con- 
ventional (OOl)-stacking of LNO/LMO SLs^- suggest- 
ing a qualitative difference between the (001)- and (111)- 
stacking directions despite having the same compositions 
and periodicity (as sketched in Fig. [IJ. This exchange 
bias in the (lll)-stacking is nontrivial since it is well 
known that LNO is a paramagnetic metal in its bulk 
formd Therefore, interesting physical questions arise: 
what is the origin of the induced magnetic moments in 
the LNO layers of these SLs? Are they in proportion 
to the charge transferred across the interfaces, consider- 
ing that electrons in the LMO layers are spin polarized? 
Why are the (001)- and (lll)-SLs qualitatively different 
with regards to their magnetic properties? And what are 
the expected results for other stacking orientations not 




FIG. 1. (Color online) The three stacking configurations of 
perovskite SLs (period n — 1) studied here: (a) (001); (b) 
(011); (c) (111). Blue and brown balls are the transition 
metals, such as Ni and Mn, while the yellow balls are the 
oxygens. The A-site cation, e.g. La, is not shown. 

yet explored, such as the (Oll)-direction? 

In this manuscript, LNO/LMO SLs will be studied 
theoretically from the perspective of microscopic models. 
Our main result is that the experimentally observed mag- 
netism in the LNO layers appears to be mainly caused by 
quantum- confinement effects. This induced magnetism is 
weakly coupled with the actual value of the charge that 
is transferred from Mn to Ni, and depends strongly and 
nonlinearly on the stacking orientations and the SL pe- 
riodicity. The underlying physical mechanism discussed 
here can partially explain the nontrivial exchange bias 
observed in (LNO)„/(LMO)„ SLs. 

Model & Method. A hybrid two-orbital tight-binding 
model is here constructed, containing the Hubbard inter- 
action for both the manganite and nickelate components 
and the double-exchange (DE) term for the manganite 
sector only. In past decades, extensive investigations 
have shown that the two-orbital DE model is a successful 
model to describe manganiteS ) 26 i 27 while the two-orbital 
Hubbard model has also been often employed for the 
nickelates . 28 ' 29 Recent theoretical studies also used these 
two models for manganite heterostructures and for LNO 
bilayersi 12 ' 30 " — For these reasons, our model provides a 
reasonable starting point to address the LNO/LMO SLs. 
More explicitly, the model Hamiltonian used here can be 
written as: 



H = H K (t ) + H™ n (J H ) + H™(V m ) + H V (U, J), (1) 
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FIG. 2. (Color online) (a) Sketch of spin-filtered quantum 
confinement. Solid and broken lines denote the spin-up and 
spin-down levels, which are degenerate for Ni but split by Jh 
for Mn. Vki < implies lower Ni levels. Pink and green 
curves are the sketch of spin-up and spin-down e g electron 
densities, (b-e) The e g -density profiles (left) and e g magne- 
tization (right) of SLs at Vm = U = 0. Black: (001)-; Red: 
(011)-; Green: (lll)-stacking SLs. Gray: LMO; White: LNO. 



where the first term is the standard nearest-neighbor 
hopping (i.e. the kinetic energy). In the results below, 
the hopping unit to is assumed to be the same for the 
bonds Mn-O-Mn, Ni-O-Ni, as well as Mn-O-Ni. This ap- 
proximation is reasonable since density functional theory 
(DFT) studies of LNO led tcM t ~ 0.6 eV while a very 
close result to ~ 0.5 — 0.6 eV was found for LMOi^i Be- 
low, to will be our energy unit. The Hund's coupling 
second term affects only the Mn ions since they are in a 
high-spin t% g configuration. Jh is here fixed to a reason- 
able value 4t (~ 2.4 eV)i2& The third term is the on-site 
e g energy of Ni, with the Mn's spin-up level as the zero of 
reference [Fig. [^a)]. This energy will be varied as a pa- 
rameter in our investigations. The last term is the stan- 
dard Hubbard interaction for multi-orbital models acting 
in the whole lattice, treated here in the Hartree mean- 
field approximation. This term contains the parameters 
U (intra-orbital Coulomb repulsion) and J (inter-orbital 
Hund exchange). The widely used ratio J = [7/4 is here 
adopted, while U will be varied as a second parameter in 
our study. More details of the model Hamiltonian and 
parameters are in the Supplementary Material. 35 

Partially due to strain, the ground state of LMO ultra- 
thin films with nearly cubic structure grown on SrTi03 
are ferromagnetic (FM)^ i 12 ' 24 In addition, e g electrons 
leak from Mn to Ni, altering the Mn valence to +(3 + 5), 
further driving the LMO layers to a FM state according 
to the phase diagram of manganitesi2& Thus, the i2 g spin 
background of the Mn layers is set here as FM unless 



explicitly noted. Then, Eq. [T] is solved self-consistently 
at zero temperature on 4 x 4 x L lattices with twisted 
boundary conditions ) 14 ' 23 where L is determined by the 
period of each SL. 

Results. First, the simplest case (U — 0, Vm = 0, n ^ 
2) is studied. This non-correlated limit (U = 0) is useful 
to clarify the underlying physics, and Vm = means that 
there is no confinement for the spin-up channel although 
the spin-down channel is still confined due to the large 
Hund Jh barrier, as sketched in Fig. [^a). The focus is 
on short periodic (n 2) SLs because it will be shown 
below that the induced magnetism is uniform in this case. 

Previous DFT calculations^ indicated that at the in- 
terface the e g electrons transfer from Mn 3+ to Ni 3+ . As 
shown in Fig. [D^b) and (d), our calculations are in agree- 
ment since even at T^i = 0, e g electrons do leak from 
Mn 3+ to Ni 3+ . This is because the e g levels in Mn's sites 
are spin polarized, pushing the Fermi level of LMO to 
higher energies since only spin-up bands can be occupied 
while spin-down bands are almost empty. By contrast, 
in LNO both the spin-up and -down bands can be filled, 
accommodating more electrons with a lower Fermi level. 

The charge transfer from LMO to LNO depends on 
both the period n as well as the stacking orientations. For 
the n = 1 case (Fig. HJ), the (001)/(011)/(lll)-stacking 
SLs have two/four/six Ni-Mn but four/two/zero Mn- 
Mn or Ni-Ni nearest-neighbors per site. For other cases 
(n 5* 2), each interfacial Mn (Ni) ion has one/two/three 
Ni (Mn) nearest-neighbors but five/four/ three Mn (Ni) 
nearest-neighbors in the (001)/(011)/(lll)-stacking SLs. 
Therefore, naively, the charge transfer from Mn to Ni 
may be the strongest (weakest) in the SLs with (111)- 
stacking [(OOl)-stacking] since they have the most (least) 
Mn-Ni bonds. However, this scenario is too simplistic. 
For example, for n — 2 the charge transferred to LNO 
is nonzero, and almost identical, for the three orienta- 
tions. However, the associated LNO magnetization is 
nearly vanishing. Moreover, in the n = 1 case the charge 
transferred to the Ni-layers is the highest for the (001)- 
stacking, yet the magnetization is the smallest for the 
same stacking. Therefore, the intensity of charge leakage 
is not in linear proportion to the induced magnetization. 

According to Fig. [21 there are several interesting fea- 
tures in the induced magnetism of LNO in these short- 
period SLs. In the n = 1 SLs, the (lll)-stacking shows 
the most prominent induced magnetization, while the 
(OOl)-stacking is nearly non-magnetic and the (011)- 
stacking interpolates between (001) and (111). In the 
n = 2 cases, the induced magnetism of Ni is almost zero 
irrespective of the stacking orientations 

These features also imply that the induced magnetism 
is indeed not simply directly correlated with the leakage 
of spin-polarized charge. For a better understanding of 
our results consider instead the band structures of our 
SLs. As shown in Fig. [3[ the spin-resolved density of 
states (DOS) at the Ni layers shows that the spin-up and 
-down channels are notoriously different. Due to the high 
spin polarization of manganites, the spin-up e g electrons 
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FIG. 3. (Color online) (a-f) Spin-resolved density of states of 
Ni. The Fermi level is at 0. Here U = and Vm = 0. Left: 
n = 1; Right: n = 2. Upper: (OOl)-stacking; Middle: (011)- 
stacking; Lower: (lll)-stacking. The e g electron densities 
and magnetism are also shown. 



spread much farther in the SLs than the spin-down e g 
electrons, which are mostly confined to the LNO layers. 
From this perspective, the LMO layers act as atomic- 
scale spin filters, causing the local band structures of the 
LNO layers to be quite different between the spin-up and 
-down channels. Furthermore, this quantum confinement 
severely depends on the stacking orientations and peri- 
odicity. The confinement is the most effective (i.e. with 
the narrowest spin-down bands) in the (lll)-stacking 
SLs due to their minimum number of Ni-Ni bonds. By 
contrast, the spin-down bands in the (OOl)-stacking are 
broader and close to the spin-up channel. 

In the n = 1 case, the narrowing of the spin-down 
bands due to quantum confinement gives rise to the ob- 
served induced magnetism, and the effect is clearly the 
most notorious for the (lll)-stacking, as discussed be- 
fore. However, since quantum confinement also exists for 
the n = 2 cases, then why their induced magnetism is 
so weak even for the (lll)-stacking? This nontrivial be- 
havior is caused by the particular quantum properties of 
this confined system. The (lll)-stacking perovskite bi- 
layer forms a honeycomb lattice that has a peculiar band 
structured In particular, for an isolated LNO bilayer 
there is a flat bottom band^Ir— In the current study, al- 
though the LNO bilayer is not isolated in the crystal, 
its spin-down channel is effectively quantum confined. 
Then, the spin-down channel has a nearly flat bottom 
band (broaden since the barrier Jh is not infinitely high) , 
that induces a large DOS peak at the band bottom [at 
—0.6 in Fig. EJf)]. This occupied localized states can 
accumulate 0.5 e g spin-down electron per Ni which sig- 
nificantly reduces the net induced magnetism. This ten- 
dency is clearly different from the n = 1 (lll)-stacking 
case, where the large DOS peak [at 0.8 in Fig. [3[c)] due 
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FIG. 4. (Color online) The Ni e g electron densities (squares) 
and magnetization (circles) versus Vsi for U — (red) and 
U = 2 (blue). Top: n = 1; bottom: n = 2. Left/middle/right: 
(001)/(011)/(lll)-stacking. 



to the confinement is actually unoccupied. In summary, 
the weak induced magnetism for the n = 2 (lll)-stacking 
is also due to the quantum confinement. 

Our results above were obtained with Vm — U = 0. 
It is necessary to confirm the robustness of those results 
with other parameters since real SLs may correspond to 
another set of V^i an d U . By varying Vjji parametric 
with U, the e g density and e g magnetization are shown 
in Fig. [4j Clearly, in all cases, the more negative Vm 
becomes, the more e g electrons accumulate on the Ni 
layers. Regarding the induced magnetism, varying Vm 
the induced magnetization remains qualitatively robust. 
In the n = 1 cases, the induced moment of Ni is quite 
weak in the (001) stacking but prominent in the (111) 
stacking, as in the (7 = case. Although the Hubbard 
repulsion U increases the induced magnetic moments, the 
qualitative features of the induced magnetization remain 
valid. Another effect of the Hubbard interaction is the 
suppression of the charge leakage from Mn to Ni. For 
example, in the n = 1 case, when U = 2 and Vm = 0, the 
e g density of Ni is close to the original value 1, imply- 
ing a very weak charge transfer. However, the induced 
magnetization remains quite prominent, further confirm- 
ing that the charge leakage from Mn to Ni is not the 
key origin of the induced magnetization in the Ni lay- 
ers. In the n — 2 case, all induced magnetic moments 
are very weak, although not exactly zero, even with the 
Hubbard interaction. All these results imply that quan- 
tum confinement effects are qualitatively robust within a 
reasonable parameter region. 

All the results above were obtained for short periodic 
SLs, in which all Ni cations are interfacial. In thicker 
cases (n ^ 3), there are inner Ni layers that do not con- 
nect directly to Mn's. Then their induced magnetism, 
if any, can be considered as a second order effect of the 
quantum confinement. As shown in Fig. the induced 
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FIG. 5. (Color online) Magnetization profiles vs. distances 
from the interfaces (e.g. the "1st" curves denote the first 
interfacial layers). A higher index denotes deeper LNO layers. 
Vm = (7 = are used. The stackings are (001) in (a), (011) 
in (b), and (111) in (c). 



magnetism in these cases displays an interesting modu- 
lation as a function of the distance from the interfaces. 37 
As expressed before, the induced magnetization in the 
(OOl)-stacking SLs is always weak, irrespective of the pe- 
riod n. Therefore, the (OOl)-stacking LNO/LMO should 
not present a robust exchange bias, in agreement with 
experiments^ By contrast, the (lll)-stacking SLs dis- 
play the largest induced magnetization, with the caveat 
that it nearly vanishes at n — 2 for the reasons already 
explained. When n ^ 3, the (111) magnetization of 
the first layer fluctuates around 0.1 ~ 0.2 /is per Ni, 
while the second Ni layers shows negative magnetization 
~ —0.1 /iB per Ni. For deeper Ni-layers (3rd, 4th, 5th, 
and more) the magnetic moments become weaker and 
weaker, finally fluctuate around zero and approaching the 
paramagnetic order of pure LNO. This qualitatively ex- 
plains the decreasing exchange bias with increasing pe- 
riod n when n ^ 5 (corresponding to the appearance 
of the 3rd layers) observed in experiments.— Although 
Fig.[5]was obtained with U = Vm = 0, other calculations 
with U — 2 and Vm = — 1 show similar tendencies with 
enhanced local moments (see Supplementary). 

Two other features of the induced magnetization are 
worth highlighting: (1) the sign oscillations with increas- 
ing distance from the interfaces and (2) the fluctuations 
of the values of the first/second layers with increasing 
period n, both suggestive of a Ruderman-Kittel-Kasuya- 
Yosida (RKKY)-likc exchange coupling between the LNO 
and LMO layers. In this sense, the almost vanishing mag- 
netism of the LNO bilayer n = 2 can also be qualitatively 
understood: the first Ni layer of the left interface is also 



the second layer counting from the right interface, lead- 
ing to a partial cancellation of the net magnetization. 

To confirm this idea, the induced magnetization was re- 
calculated by flipping the magnetic background of some 
Mn layers in the (lll)-stacking n — 2 case. Two situa- 
tions were tested: (1) flipping one layer in each bilayer 
(Mnt-Mn|-Ni-Ni-Mnt-Mn|); (2) flipping one bilayer en- 
tirely (Mn^-Mn|-Ni-Ni-Mnt-Mnt). Both these cases now 
give very prominent induced local magnetization in the 
LNO n = 2 bilayers, which couples antiferromagnetically 
between these two neighboring layers (-NiJ.-Nit-). While 
this result is different from the ferromagnetically ordered 
case (Mnt-Mnt-Ni-Ni-Mnf-Mnf) case, it also supports 
the scenario of a RKKY coupling between the Mn's local 
moments and the Ni's induced moments. 

It is interesting to compare the induced magnetism of 
the LMO/LNO SLs against the magnetic reconstruction 
observed in the LMO/SMO SLs that has been intensively 
investigated £~— The induced magnetism in pure man- 
ganites SLs is associated with the magnetic phase tran- 
sitions driven by the modification in the e g density and 
straini 12 i 14 ' 38 Although not precisely in a one-to-one cor- 
respondence, in pure manganite SLs the induced mag- 
netism can be traced back to the phase diagram of bulk 
manganites. By contrast, in the LMO/LNO SLs, the 
induced magnetization, which is nearly decoupled from 
the e g electron transfer, is qualitatively different from 
the pure manganites SLs, but is more analogous to the 
dimensionality control of electronic phase transitions in 
LNO-LaA10 3 heterostructures3^ 

Conclusions. The induced magnetization found in 
LNO in the LNO/LMO SLs was studied here via a hy- 
brid microscopic model. The results of our model agree 
with previous DFT investigations^ but provide addi- 
tional details and a deeper physical insight. Summariz- 
ing our conclusions: (1) in the n — 2 SLs with FM LMO, 
the LNO layers are nearly non-magnetic independently 
of the stacking directions. (2) The induced magnetiza- 
tion of LNO in the (lll)-stacking SLs is always the most 
prominent. By contrast, the LNO layers in the (001)- 
stacking SLs are always nearly non-magnetic, compati- 
ble with the exchange bias investigations.— The results 
for the (Oll)-stacking are in between those of the (111)- 
and (OOl)-stackings which can be verified in future exper- 
iments. (3) The induced magnetic moments of the first 
Ni layers are parallel to the moments of their nearest- 
neighbor Mn layers, but the second Ni layers usually dis- 
play negative moments if not zero. The induced local 
magnetism decreases to zero in an oscillatory manner by 
increasing the thickness of the LNO layers. The under- 
lying physical mechanism for the induced magnetization 
is associated with the spin-filtered quantum confinement 
supplemented by a RKKY-like exchange coupling, qual- 
itatively different from the magnetization reconstruction 
in most previously studied oxide heterostructures. 
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